ON A QCD-BASED PION DISTRIBUTION AMPLITUDE VS. RECENT 

EXPERIMENTAL DATA 



A.P. BAKULEV 1 ' a , S.V. MIKHAILOV 1 , N.G. STEFANIS 2 
1 Bogoliubov Laboratory of Theoretical Physics, JINR, 141980 Dubna, Russia 
2 Institut fur Theoretische Physik II, Ruhr-Universitdt Bochum, B-44180 Bochum, Germany 




Using QCD sum rules with nonlocal condensates the twist-2 pion distribution amplitude is 
determined by means of its moments and their confidence intervals, including also radiative 
corrections. An admissible set of pion distribution amplitudes is constructed in the £12, 04 plane 
of the Gegenbauer polynomial expansion coefficients. The determined 02, 14 region strongly 
overlaps with that extracted from the CLEO data by Schmedding and Yakovlev. Comparisons 
are given with results from Fermilab experiment E791 and recent lattice calculations. 



1 Pion Distribution Amplitude 

The main object of this talk is the pion distribution amplitude (DA), defined by 

(0 I d(z) lul5 E(z, 0)u(0) I tt(P)} = iUP u f dx e tx ^ <p n (x, M 2 ) , (1) 

* -0 Jo 

where gauge invariance is ensured by the Fock-Schwinger string E(z, 0) = V exp \ig J Q A^tjdt^ . 

The pion DA has the following important properties: (1) it is multiplicatively renormalizable, 
(2) it has isospin symmetry: (f n (l — x,/j?) = ip n (x,^i 2 ), and (3) its normalization is conserved: 
Jq 1 dx if n (x,^ 2 ) = 1. 

Pion DAs naturally appear in the perturbative QCD (pQCD) description of any hard exclu- 
sive process with pions. For example, the form factor of 7*7* — > vr° decay with — qf 2 ~ Q 2 > 
1 GeV 2 is factorized in pQCD according to 

F rr ^ n o(qlqi) = C{qlql, f i 2 ;x)®ip n (x; f i 2 )+0(Q-*) +... (2) 
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For (p n (x) we use an expansion in terms of eigenfunctions of the 1-loop evolution kernel, 

XxCn /2 (2x- 1), 

^(z; = <^ s (x) \l + a 2 {f)Cl /2 {t) + a 4 (^)Ct /2 (0 + . . 



(3) 



with </^ s (x 



6xx being the asymptotic pion DA and £ = 2x — 1. In this expansion all scale- 
dependence is accumulated in the coefficients {ao(u 2 ), a 4 (fj, 2 ), . . .}. Note that the evolution of 
the pion DA at the 2-loop level is available from Era 

How one can obtain the <£v(x, /x 2 )? It is possible to extract it from: 

• experimental data: (i) see the recent papers of the CLEO Collaboration! and the analysis of 
Schmedding and Yakovlev of these dataQ, and (ii) using the data of the E791 Collaboration! 

• QCD Sum Rules with Non-Local Condensates (NLC) see BOA 

• transverse lattice simulations!!^ 

• instanton-induced models EM 2 ]. 

In this talk we consider all these sources separately but the main focus is on the first 2 items. 



2 Revision of the NLC QCD SR Results 



We re-analyze our NLC SRs with the modification of one of the antiquark-gluon-quark NLC 
contributions to obtain revised values of the moments (^ Ar ) 7r = J <Ptt(x; fi 2 ) (2x — 1) N dx (N = 
2, . . . , 10), new estimates of error-bars, and a new estimate of (x~ 1 )sr = Jq ip n (x; fi 2 )x~ 1 dx, 
where a SR is used that is constructed directly for this quantity. 

Our model of NLCs, illustrated by the scalar NLC (q(0)q(x)) = (g(O)g(O)) exp (-|x 2 |A 2 /8), 
uses only a single parameter A 2 , which is related to the vacuum (fields) correlation length 
1/Aq = A the latter being of the order of the hadron size, as estimated in non-pQCD approaches: 

0.4 ± 0.1 GeV 2 [ QCD SRs, 1983 i£ 

0.5 ± 0.05 GeV 2 [ QCD SRs, 1991 J^L (4) 
« 0.5 GeV 2 [ Lattice, 1998-99 B@ ] 

From the lhs of Fig. 1, one sees that the quality of the stability in the NLC QCD SRs for 
2 = 0.4 GeV 2 is quite high (solid line stands for the best threshold sq = 2.2 GeV 2 , dashed lines 
for 10%-variations of this parameter). The obtained moments are shown on the lhs of Fig. 2. 
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Figure 1: Leftside: stability curves for {a; -1 )^ against the Borel parameter M 2 . Central part: Bunch of admissible 
DAs corresponding to X 2 , = 0.4 GeV 2 with best-fit parameters ai = +0.188, = —0.130. Right side: Same as 
central part, but for A 2 , = 0.5 GeV 2 with best-fit parameters a2 = +0.126, 04 = —0.091. 

We reconstruct the pion DA from these five moments, using models (|3[) at /x 2 = 1 GeV 2 with 
two non-zero Gegenbauer coefficients. The best-fit DAs obtained this way (with x 2 ~ 10~ 3 ) 
are shown in Fig. 1 as thick solid lines. The corresponding error bars to the DAs, allowed 
by the moment SRs, are also shown. The bunches of these broken lines represent the self- 
consistent DAs in the sense that the value of the associated inverse moment, (x" 1 )^ = 3.17(8) 
((a^^Tr = 3.13(8), for A 2 = 0.5 GeV 2 ), is in good agreement with the value determined from 
the special SR: (x^ 1 )^ = 3.33(32) ((x" 1 )^ = 3.19(32)). 



3 Comparing with CLEO results in SY approach 

Schmeddineand Yakovlevi have provided a useful analysis of the CLEO data on the 7* (5)7 — > tt 
form factor H, using light-cone QCD SRs and taking into account NLO and twist-4 corrections. 
They estimated the first two Gegenbauer coefficients of the pion DA and obtained 



a 2 = 0.19 ± 0.04(stat) ± 0.09(sys) , a 4 = -0.14 ± 0.03(stat) ^ 0.09(sys) 



(5) 



with results displayed in Fig. 2 in the form of confidence regions in terms of a2 and 04. 

We evolve our allowed sets to the CLEO scale // 2 = (2.4 GeV) 2 and insert them directly 
into the SY plot to obtain results in the a2,&4 plane, shown in Fig. 2. Inspection of these plots 
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Figure 2: Left: (£ N ) against TV. Central part: Confidence region of pion DAs for A 2 = 0.4 GeV 2 . Right side: the 
same as before with specific explanations given in the text, but with A 2 = 0.5 GeV 2 . 

reveals that bunch-1 [A 2 = 0.4 GeV 2 ] and bunch-2 [A 2 = 0.5 GeV 2 ] are intersecting with the 
SY lex-region and are thus quite compatible with CLEO data. 

Inspecting the SY plot, one natural question arises: why the confidence regions are stretched 
along the diagonal a 2 + 04 = const ? To answer this question, it is instructive to analyze the 
CLEO data, employing pure pQCD: 



3Q^ 

47T 



Fr 



77 7T 



Q 2 -C®(p n = (x~ 1 } 7r + 3Q s (Ao + A2 + A 4 )-r- 



(6) 



We see, that, up to radiative corrections, CLEO in fact measured the inverse moment of the pion 
DA, (x -1 )^, which is simply connected to the diagonal combination of the a 2 and 04 coefficients: 



3 (1 + a 2 + a 4 ) 



(7) 



The SY analysis gives: a 2 + a 4 = 0.05 ± 0.07, whereas the special NLC QCD SRs yields for this 
moment i(x _1 )^ R - 1 = 0.10 ± 0.10 (at y. = 1 GeV), and our bunch-1 of the allowed pion DAs 
produces a 2 + a 4 = 0.056 ± 0.03 (at \i = 2.4 GeV), in excellent agreement with S&Y result. 
It is useful to have a look on the numerical values of different terms in (S) 



3Q 2 

47T 
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3 [1 + (02 + a 4 ) + a s A + a s (A 2 + A 4 ) + . . . ] 
3[1+ 0.05 - 0.17 - 0.014 + ... ] 



(8) 



The CLEO data gives numerically (3Q 2 /4tt) F^* n « 2.45. We see that in the LO pQCD analysis 
(i.e., when all a s A]\f = 0), one arrives at the estimate: (x~ 1 ) 7T = 2.45. From this, one might 
conclude that (p n (x) should be narrower than the asymptotic one0. Taking into account only 
the main part of the NLO correction (a s Ao = —0.17), we get the estimate, (x -1 )^ = 2.96, 
and conclude that ip n (x) could have the same width as the asymptotic DA. But the full NLO 
(plus twist-4 contribution) light-cone QCD SR§ provides instead, (x 1 } 7r ~ 3.15, indicating that 
(p n (x) should be broader than the asymptotic DAt3, just as it appears in our bunches in Fig. 1. 

The E791 collaboration has measured diiet production in diffractive ttA interactions!!. Such 
an experiment has been suggested in 1993 inll^l as a means of measuring directly the squared pion 



DA at large transverse momentum transfers. We obtain a good fit of the E791 data using our 
model (symbol * on the rhs of Fig. 2) with = 0.12 and af* = 0.01 at the scale ~ 8 GeV 2 . The 
resulting value of the "diagonal" at the CLEO scale appears to be too large: + af* — 0.14. 
In our view, the interpretation of this experiment is still questionable. Moreover, it seems that 
errors are too large and should be estimated more carefully. 

There are two recent papers involving transverse lattice simulations. Dalleyi produced (see 
symbol A on the rhs of Fig. 2) 

^{x-y 2 ~ 1 GeV 2 ) = 6xz[l + 0.133<?2 /2 (2x - 1)1 , (9) 

and, on the other hand, Burkardt and Seal0 arrived - using the same approach - at a different 
DA (denoted by • in Fig. 2) very close to the asymptotic pion DA. Note that this large difference 
seems to indicate that the errors of this method are still large and should be estimated more 
precisely. 

The existing predictions from instanton-induced models EJO are too close to ip^, except for 
the new model by Praszalowicz and RostworowskiE3, which is just outside the confidence region 
of the (/32-bunch and on the boundary of the 95%-region of SY (symbol ♦ on the rhs of Fig. 2). 

Acknowledgments 

This work was supported in part by the Russian Foundation for Fundamental Research (contract 
00-02-16696), INTAS-CALL 2000 N 587, the Heisenberg-Landau Program (grants 2000-15 and 
2001-11), and the COSY Forschungsprojekt Julich/Bochum. 

References 

1. S. V. Mikhailov and A. V. Radyushkin, Nucl. Phys. B 273, 297 (1986). 

2. D. Midler, Phys. Rev. D 49, 2525 (1994); Phys. Rev. D 51, 3855 (1995). 

3. J. Gronberg et al. (CLEO Collaboration), Phys. Rev. D 57, 33 (1998). 

4. A. Schmedding and O. Yakovlev, Phys. Rev. D 62, 116002 (2000). 

5. E. M. Aitala et al. (Fermilab E791 Collaboration), |hep-ex/0010~043l ; 

D. Ashery, |hep-ex/9910~024] ; Nucl. Phys. Proc. Suppl. B 90, 67 (2000). 



6. S. V. Mikhailov and A. V. Radyushkin, JETP Lett. 43, 712 (1986). Sov. J. Nucl. Phys. 
49, 494 (1989); Phys. Rev. D 45, 1754 (1992). 

7. A. P. Bakulev and S. V. Mikhailov, Z. Phys . C 68, 451 (199 5); Mod. Phys. Lett. A 11, 
1611 (1996); Phys. Lett. B 436, 351 (1998); jhep-ph/0006206| . 



8. A. P. Bakulev, S. V. Mikhailov and N. G. Stefanis, |hep-ph/0103lT9 



9. S. Dalley, |hep-ph/0101318 



10. M. Burkardt and S. Seal, |hep-ph/010l338 



11. I. V. Anikin, A. E. Dorokhov, and L. Tomio, Phys. Part. Nucl. 31, 509 (2000). 

12. V. Yu. Petrov et al., Phys. Rev. D 59, 114018 (1999). 

13. V. M. Belyaev and B. L. Ioffe, Sov. Phys. JETP 57, 716 (1983); 

A. A. Ovchinnikov and A. A. Pivovarov, Sov. J. Nucl. Phys. 48, 721 (1988). 

14. A. A. Pivovarov, Bull. Lebedev Phys. Inst. 5, 1 (1991). 

15. M. D'Elia, A. Di Giacomo, and E. Meggiolaro, Phys. Rev. D 59, 054503 (1999). 

16. E. Meggiolaro, Nucl. Phys. Proc. Suppl. 83, 512 (2000). 

17. A. V. Radyushkin and R. T. Ruskov, Nucl. Phys. B 481, 625 (1996). 

18. N. G. Stefanis, W. Schroers, and H.-Ch. Kim, Eur. Phys. J. C 18, 137 (2000). 

19. L. Frankfurt, G. A. Miller, and M. Strikman, Phys. Lett. B 304, 1 (1993). 

20. M. Praszalowicz and A. Rostworowski, Cracow preprint , TPJU-4/2001 (2001). 



